EXPLICIT UPPER AND LOWER BOUNDS 
FOR THE TRAVELING WAVE SOLUTIONS 
OF FISHER-KOLMOGOROV TYPE EQUATIONS 



A. GASULL*, H. GIACOMINI, AND J. TORREGROSA* 



Abstract. It is well-known that the existence of traveling wave solutions for 
reaction-diffusion partial differential equations can be proved by showing the 
existence of certain heteroclinic orbits for related autonomous planar differential 
equations. We introduce a method for finding explicit upper and lower bounds 
of these heteroclinic orbits. In particular, for the classical Fisher-Kolmogorov 
equation we give rational upper and lower bounds which allow to locate these 
solutions analytically and with very high accuracy. 



I. Introduction and Main Results 
Consider the adimensionalized reaction-diffusion partial differential equation 

u t = u xx + f(u), (1) 

of Fisher-Kolmogorov type, where f(u) is a smooth function satisfying certain 
hypotheses. 

The usual Fisher-Kolmogorov equation corresponds to f(u) = u(l — u) and mod- 
els the spreading of biological populations, see [21 [5] . Other well-known cases are 
the Newell- Whitehead- Segel equation, f(u) = u(l — u 2 ), for describing Rayleigh- 
Benard convection, see [HI E|, and the Zeldovich equation, f{u) = u{l — u){u — a) 
with < a < 1, that appears in combustion theory, see [IT]. See also [H U\ I5| HP]. 

It is known that the traveling wave solutions u = u(x — ct) of satisfying 

lim u{s) = and lim u(s) = 1, 

s—t—ao s— >oo 

appear when u(s) is a special solution of the second order equation 

u + cu + f{u) = 0, 

where the dot indicates the derivative with respect to s. This solution can be seen 
as the heteroclinic orbit H c of the planar system 

■ = W ' f( 1 ( 2 ) 
v = —cv — jyu), y ' 

that connects the origin with the saddle point at (u,v) = (1,0), see Figured] For 
instance for the classical Fisher-Kolmogorov equation it exists only when c > 2. 
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Figure 1. Phase portrait showing the unstable separatrix H c . 



The goal of this paper is to give analytic upper and lower bounds of the het- 
eroclinic orbit H c as well as of their time parametrization. We will approach to 
this question with similar tools to those introduced in [3] . A key point consists in 
using the local behavior of the separatrices of the critical points to guess global 
algebraic bounds for the actual orbits. 

First we prove a general result for system (|2J) . It is illustrated in Figure |2J 

Theorem 1.1. Consider system (J2J), with f satisfying /(0) = /(l) = 0, /'(0) > 0, 
f'{l) < 0, f"(u) < for allu G M and c > 2y/f'(0). Let H c be its heteroclimc 
orbit and define 



Then H c can be parametrized as H c = {(u, h c (u)), u e [0, 1]} and for allu £ (0, 1), 

-A/(u)</i c (m)<-A/(m)<0. 

Moreover, if (u c (s),v c (s)) is the parametrization of H c such that u c (0) = 1/2 it 
holds that 



where (a, b) denotes the smallest closed interval containing a and b and u = z\(s) 
is the solution of the Cauchy problem 




and X : 



c — 



2/'(l) 



u c (s) £ (zx(s),zj(s)) , 




V 



u 




v = -\f(u) 



u 




s 



(i) Phase pi; 



anc 



il) Time parametrization 



Figure 2. (i) Upper and lower bounds of the unstable separatrix 
H c (dotted line), (ii) Bounds for the time parametrization of H c . 
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For instance when we particularize the above theorem to the classical Fisher- 
Kolmogorov system we obtain that 

u c (s) £ / — -^V (3) 
V ; \ 1 + eh s 1 + e Xs I 

The results obtained for system (j2J) can be improved when we study this case. We 
prove: 

Theorem 1.2. Let H c = {(u, h c (u)), u G [0, 1]} be a parametrization of the hete- 
roclinic solution of system (j2J) when f(u) = u(l — u). For u G (0, 1) and c > 2 it 
holds that 

h c 2 (u) < hl(u) < ■■■ < h c 100 (u) < h c {u) < R c 10 (u) < Rg(u) < ■ ■ < R\{u), 

where h^{u) is the Taylor polynomial of degree n of the unstable separatrix of 
the saddle point (1,0) of the system and each R c m is a rational function whose 
numerator has degree m + 2 and its denominator degree m, constructed from the 
Pade approximants of v = h2 m+2 (u). See Figure^ 

* v 




Figure 3. Upper and lower bounds of the separatrix H c of the 
Fisher-Kolmogorov equation, which is plotted as a dotted line. 

For a given value of c it is also possible to study the maximum distance be- 
tween two of the above functions. For instance when c = 99/100, we can prove 
that -Rg 9 ^ 100 (x) — /i2(f 10 °( x ) < 2 x 10~ 19 . It is also possible to go further in the 
computations. For example we get that 

C 10 » < />" /10 V) < < A >) and i4 9/1 °°(«) - hT°°(u) < 2 x 10" 47 , 

for all u G [0, 1] and the maximum error is at u = 0. For bigger values of c it is 
needed to compute approximations of higher degree to arrive to similar bounds of 
the error. 

The s-parametrization of H c can also be obtained with more accuracy for the 
classical Fisher-Kolmogorov case. We only present here a first result. Sharper 
approximations are detailed in Section HI 

Theorem 1.3. Let (u c (s),v c (s)) be the time-parametrization of the heteroclinic 
orbit H c of system ([2]) with f(u) = «(1 — u) and such that u c (0) = 1/2. Define 

3 + 2V2 



U' 



2 ■ 



(1 + v ^+ e (v^+4-^/2)' 
Then: 

(i) When 2 < c < 5/y/E it holds that sgn(w c (s) — U c (s)) = — sgn(s) 
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(ii) When c = 5/ a/6 it holds that u c (s) = U c (s). 
(Hi) When c > 5/\/6 it holds that sgn(u c (s) — U c (s)) = sgn(s). 

The above inequalities improve the bounds given in OH]). Notice also that when 
c = 5/\/6 we have obtained the exact expression of the traveling wave solution of 
the classical Fisher-Kolmogorov equation, u t = u xx + u{u — 1), 

3 + 2y/2 

U{x,t) = — — — , 

1 + \/2 + e7SV x "75*J 

which coincides with the one given in [T] . The novelty of our result is that a similar 
expression gives a bound of the traveling wave for all the values of c. 

The methods developed to study the classical Fisher-Kolmogorov equation can 
also be applied for the Newell- Whitehead- Segel and the Zeldovich equations. 

2. Proof of Theorem 11.11 

For computational reasons it is more convenient to locate the saddle point of 
system ([2]) at the origin. So we introduce the new variables x = 1 — u, y = v and 
t = s and it writes as 

X ' = ~ Vj U) 

y = -cy + g(x), 

where g(x) = — /(l — x) and the prime denotes derivative with respect to t. 
Observe that the variable x introduced above does not coincide with the one used 
in equation ([T]). Notice that g satisfies the following set of hypotheses 
H : 0(0) = 0(1) = 0,g'(0) < 0,g'(l) > and g"(x) > for all x G R. 

The above system has only two critical points (0,0) and (1,0). The origin is a 
saddle point with eigenvalues 



AJ = - C± ^- 4g ' (0> (5) 

and corresponding eigenvectors (1, — A+) and (1, — Aj). Similarly, when c 2 — 
4</(l) — tli e point (1,0) is an attracting node and its eigenvalues are 



2 

A sketch of the phase portrait of system (j3J) is given in Figure HI There we can 
see the heteroclinic connection that we are interested to locate. Indeed it is given 
by ones of the branches of the unstable separatrix of the saddle point. In the new 
coordinates we will call it r c . 

Proposition 2.1. Consider system (j3J), with g satisfying hypotheses H and c > 
2yM)). Set 



C- y/c2-4g'(l) , T C- y/c 2 -4^(0) 

A := — and A := — . 

2g'(l) 2g'(0) 

LetT c be its heteroclinic orbit. Then T c can be parametrized as T c = {(x, ^y c (x)), x G 
[0, 1]} and it holds for all x G (0, 1) that 

Xg(x) < 7 c (x) < Xg(x) < 0. 
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Figure 4. Phase portrait showing the unstable separatrix T c of system (j2J). 
Proof. Consider the 1-parameter family of maps G\ (x, y) = y—Xg(x). We compute 



(VG A (x, y), ( - y, -cy + g(x)))\ = g(x) (l - cA + g'(x)X 2 ) =: g(x)N x (x). 



Since on (0, 1) it holds that g(x) < 0, if we choose A such that N\ does not 
vanish on the same interval we will have that the corresponding curve y = Xg(x) 
is without contact. Notice that N' x (x) = X 2 g"(x) and so for A ^ 0, the hypotheses 
H imply that the function N\ is increasing. Therefore: 

(a) If for some A > 0, iV A (0) > then N x (x) > for all x > 0. 

(b) If for some A > 0, iV A (i) < then N x (x) < for all x < 1. 

The conditions N x (j) = for j = 0, 1, write as 1 — cA + g'(j)\ 2 = 0. Their 
solutions are 



Therefore, taking Gj(x, y) = and G\(x, y) = as an upper and lower bound- 
ary, respectively, we have constructed a subset of the strip {(x,y), < x < 1}, 
that contains the heteroclinic orbit T c , see Figure El Then, on it, x' = —y > 0, T c 
can be parametrized as a function of x, say y = 7 c (x), and the inequalities of the 



C± y/c 2 -4^(j) 

2g'(j) 



It is easy to prove that 



A+ < < A, 



: A < A 




statement follow. 



□ 



y 



X 




Figure 5. Upper and lower bounds of the unstable separatrix (dot 
ted line) Y c . 



Notice that y = Xg(x) = — X^x + 0(x 2 ) and therefore this curve is tangent to 
the unstable separatrix of the saddle point. 
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Using the above proposition we can also approach the parametrization of r c 
with respect the actual time t. We need to introduce some new functions. Given 
the Cauchy problem 

t = - X9{x) > (6) 

x(0) = 1/2, 

we will denote by x = w\(t) its corresponding solution. 

Theorem 2.2. Under the same hypotheses and notations of Proposition \2.1\ if 
{x c (t),y c (t)) is the parametrization ofT c such that x c (0) = 1/2 it holds that 

x c {t) e (wx(t),wj(t)) , 

where W\(t) is defined in (JB]) ; see Figure® 




Figure 6. Upper and lower bounds of the t-parametrization of T c 
(dotted curve). 

Proof. By Proposition 12.11 we know that 

Xg(x) < 7 c (x) < Xg(x) < 0. 

Recall that T c is parametrized as y — j c (x) and therefore since in system (jl]), 
x'(t) = —y(t), it holds that x' c (t) = — 7 c (x c (t)). Therefore 

-Xg(x c (t)) > x' c (t) > —Xg(x c (t)) > 0. 

Since x c (0) = wj(0) = 1/2, x c (t) satisfies the differential inequality 

x' c (t) > -Xg(x c (t)) 

and wj(t) the equality, it holds that x c (t) > wj(t) for all t > 0, as we wanted to 
see. The other cases follow similarly. □ 

Corollary 2.3. // in system (jl]) , g(x) = x m — x, 2 < m G N, then 

1 1 

x c (t) G 



In particular for the Fisher-Kolmogorov case, m = 2, 
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Proof. It suffices to solve the Cauchy problem ([6]) and then apply Theorem 12.21 
When g{x) = x m — x we obtain 

wx(t) = 



_|_ - l) e -(m-l)Xt 

and so the result follows. □ 

Theorem 11.11 is a reformulation for system (J2J) of the results of this section 
obtained for system (j3J). 

3. Sharper upper and lower bounds for r c in the 

FlSHER-KOLMOGOROV CASE 

This section will be devoted to find sharper upper and lower bounds for T c in the 
Fisher-Kolmogov system when g(x) = x(x — 1). We will use dynamical tools 
inspired in [3J. One of the key points will be to find algebraic curves constructed 
by imposing that these curves coincide as much as possible with the unstable 
separatrix of the saddle point. 

To avoid the appearance of square roots during the computations it is convenient 
to include a new parameter r in such a way that 



1 

r. 

r 



Then system (jl]) writes as 



x = -y, 

1 



y' = (r )y + x(x-i), re(0,\/2-l]. 



r 

Notice that the condition on r implies that c = r — 1/r > 2. One advantage 
of introducing this new parameter is that the eigenvalues of the saddle are now 
— 1/r < < r. In the notation of the previous section A+ = A = r. Hence, if we 
denote T r := Ti/ r _ r the searched heteroclinic trajectory, from Proposition 12 . 1 1 we 
know that y = rx(x — 1) is an upper bound for T r . 

First we need to know the local expansion of the unstable manifold of the saddle 
point. 

Lemma 3.1. The local unstable manifold of the origin of system ([5D writes as the 
analytic function 

y = n [x)= — rx-\ x + -. — - — -^-r-. — rx 

y K J 2r 2 + 1 (2r 2 + l) 2 (3r 2 + 1) 

10r 5 4 

+ (2r 2 + l) 3 (3r 2 + l)(4r 2 + if 

12r 7 (19r 2 + 6) 5 5 



(2r 2 + l) 4 (3r 2 + l) 2 (4r 2 + l)(5r 2 + 1) 
and the subsequent terms can be determined recurrently and are positive for r > 0. 
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Proof. Let y = h r (x) = h(x) be the local expression of any of the separatrices of the 
saddle point, being h an analytic function at zero. Then y' — h'(x)x' 
or equivalently, 



y=h(x) 



= o, 



r ) y + x ( x — 1) + h'{x)y 



y=h(x) 



r - - + h'{x) h(x) - x + x 2 = 0. (9) 



Writing = h\X + /i 2 x 2 + . . . and plugging this expression in the above one we 
get that h\ + (r — l/r)hi — 1 = 0. So we choose /ii = — r. Then the right hand 
identity in Q writes as 



h 2h 2 x + 3/130^ 

r 



-rx + /i 2 x 2 + h 3 x 3 + 



x + x =0. 



So — 2/i 2 rx 2 — h 2 x 2 /r + x 2 = 0, which implies h 2 = r/(2r 2 + 1). In general, for 
n > 2 it holds that 



— /i„ + 2h 2 h n ^ 1 + 3h 3 h n ^ 2 H + (n — l)/i n _i/i 2 - ft./^ = 0. 

r 



Therefore, 



2/i 2 /i n -i n—2 



(n - l)h n -ih 2 



nr + 1 /r 

and by induction /i„ > 0, for > 2, as we wanted to prove. 



□ 




Figure 7. Upper and lower bounds of the separatrix T c plotted as 
a dotted thin line. Here 2 < n < 100 and 1 < m < 10. 

Proposition 3.2. Lei T r = {(x,7 r (x)), x G [0,1]}, 6e a parametrization of the 
heteroclinic solution of system Then for all x G (0, 1) it holds that 

h r 2 (x) <h r 3 (x)<---< h r 99 (x) < h r 100 (x) < 7 r (x), 

where h r n (x) is the Taylor polynomial of degree n in powers of x of the function 
h r (x) defined in Lemma \3.1\ See Figure^ 

Proof. We will give first the details for n = 3. It is evident that for x > 
sufficiently small we have h r (x) > h r 3 (x). First we evaluate the polynomial h 3 (x) 
at x = 1. We obtain 

2r 5 (5 + 6r 2 ) 
31 ' (l + 2r 2 ) 2 (l + 3r 2 )' 
This quantity is negative for r > and it is a monotonous decreasing function 
of r. Its value at r = a/2 — 1 is approximately —0.054. Moreover, /^(O) = 0, 
h 3 (x)" > for x > and ^3(2) < for x > and sufficiently small. Then we 
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conclude that h r z {x) < for all x G (0, 1]. We will show now that the flow of the 
vector field Z r (x, y) = (—y, (r — l/r)y + x(x — 1)) associated to system (jSJ) crosses 
the curve F£(x, y) = y — h r 3 (x) = upwards. We compute 

M;{x) := (VFax,y),Z r (x,y))\ y=h r (x) . 

We obtain 

10rV 12r 6 x 5 



m;(x) 



'1 + 2r 2 ) 3 (l + 3r 2 ) (1 + 2r 2 ) 4 (l + 3r 2 ) 2 ' 



This expression is positive for x > 0. 

All the other cases can be studied by using the same method. The key point is 
that all the monomials in the corresponding expression of M*(x) are positive for 
x > 0. □ 



In Proposition 13.21 we have chosen to stop at the value 100 because the compu- 
tation of the function M[ 00 used in its proof together with the testing that all its 
coefficients are positive takes more than four hours of CPU time in our computer. 
In any case, for practical uses, it suffices to consider small values of n. 

Remark 3.3. It is clear that if we could prove that the radius of convergence of 
the series given in Lemma \3.1\ when r G (0, v2 — 1], is 1 we would have obtained 
an infinite monotonous sequence of polynomials tending to the actual separatrix. 
Unfortunately we have not succeeded in our attempts. In any case, computing an 
approximation of the radius of convergence, by using several hundreds of terms of 
the series, seems to show that the result is true. 

We give now upper bounds for r r . Before to state the next proposition we 
introduce some definitions. Consider the rational function 

R r (x)- ^) 



rx(x — 1) 

We define the sequence of rational functions R r n {x) as follows: 

R r n (x) = rx(x - l)P; n (x), n = 1, . . . , 10, 

where the P^ n {x) are the Pade approximants of order (n, n) of the function R r (x). 
For instance 

(2r 2 + l)(3r 2 + l)-3r 2 a; 
l{X> rX{X j (2r 2 + l)(3r 2 + l)-r 2 (6r 2 + 5)a;' 

Recall that the Pade approximants P£ n (x) are rational functions whose numera- 
tors and denominators are polynomials of degree n and their Taylor expansions in 
powers of x are the same that the Taylor expansion of R r (x) up to order 2n. We 
write 

TO= ra (x-l)^|, with C7;(0)>0. 

The rational functions R r n (x) vanish at x = and x = 1 and their Taylor expan- 
sions coincide with the Taylor expansion of h r (x) up to order In. As we will see 
in the proof of next proposition they are well defined for all x G [0, 1] because 
CUx) ± 0. 
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Proposition 3.4. Let T r = {(x, y(x)), x G [0,1]}, be a parametrization of the 
heteroclinic solution of system (jHJ). For x G (0, 1) and r G (0, a/2 — 1] it holds that 

7 r (x) < R r 10 (x) < R r 9 (x) < ■ ■ ■ < R r 2 (x) < R[(x), 

see Figure^ 



Proof. We will show that for each n, 1 < n < 10, the flow of the vector field 
Z r (x, y) = (—y, (r — l/r)y + x(x — 1)) associated to system (|S]) crosses each curve 
G r n {x,y) = y — R n (x) = forwards. We compute 



(VG r n (x,y),Z r (x,y))\ y=Rrn{x) 



We obtain 



n: 



x 



r 4n+2 a ,2n+2 ( 



b:(x) 



X 



where B r n (x) and C^(x) are polynomials of degree n in x. The coefficients of B r n (x) 
and C^(x) are polynomials in r. 

We want to prove that for r > both polynomials B n (x) and C n (x) are positive 
for x G (0, 1). We will approach the problem varying the parameter r and studying 
which are the possible bifurcations for the number of zeros of these polynomials 
when x G (0, 1). The values of B n (0), C n (0), B r n {l) and C£(l), are polynomials 
in r whose coefficients are natural numbers. Then we have -B^(O) > 0, C^(0) > 0, 
5^(1) > and C^(l) > for r > 0. For n > 2, we evaluate now the discriminants 
of the polynomials B^(x) and C^(x), Dis(5^(x), x) and Dis(C^(x), x). We obtain 



a polynomial in r of degree b n with positive integers numbers as coefficients for 
Dis(5^(x), x) and a polynomial in r of degree c n with also positive integer numbers 
as coefficients for Dis(C^(x), x), where the degrees b n and c n are given in Table El 
Then, the two discriminants do not vanish for r > 0. In consequence, the number 
of real roots in (0,1) of each polynomial, B r n (x) or C^(x), does not change for 
r > 0. Picking a concrete value of r, for instance r = 1/10, we found, by applying 
the Sturm algorithm, that they have no roots in (0,1). Then, we deduce that 
B^x) > and C^(x) > for x G (0, 1) and r > 0. Therefore, we have proved 
that for n > 2, N^(x) < for x G (0, 1) and r > 0, as we wanted to see. The case 
n — 1 is much easier because we do not need to compute the discriminants and 
we omit the details. 



n 


2 


3 


4 


5 


6 


7 


8 


9 


10 


b n 


40 


212 


624 


1480 


2900 


5028 


8260 


12560 


18180 




24 


100 


264 


584 


1100 


1860 


2996 


4496 


6444 



Table 1. Degrees of the numerator and the denominator of the 
resultants computed in the proof of Proposition 13.41 



Finally for k = 2, 3, . . . , 10 we obtain by a direct computation 

x 2k {x- l)r 4k - l D k {r 2 ) 



R r k (x)-RU(x) 



Cl^Ci^x) 



where Dk(x) is a polynomial with positive coefficients. Hence R r k {x) < i?^_ x (x) as 
we wanted to show. □ 
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As a corollary of Propositions 13.21 and 13.41 we obtain the following result: 

Theorem 3.5. Let T r = {(x, 7 r (x)), x G [0, 1]}, be a parametrization of the hete- 
roclinic solution of system (jSJ). For x G (0, 1) and r G (0, y/2 — 1] it holds that 

h r 2 (x) < hl(x) < ■■ ■ < h r 100 (x) < 7 r (x) < R[ (x) < R r g (x) < ■■■ < R\(x). 

If we were interested in obtaining more precise upper bounds for the heteroclinic 
orbit, we could simply increase n and apply the same procedure. For instance for 
r = 1/10 we have performed all the computations and proved that for x G (0, 1), 

h% w (x) < 7 1/10 (x) < Rli W (x) and R& 10 (x) - h\{ 10 (x) < 9 x 1(T 32 , 

hT(x) < 7 1 / 10 (x) < rT(x) and - ttf°(x) < 2 x 10" 40 , 

hl{ 10 (x) < 7 1 / 10 (x) < i?^ 10 (x) and R% W {x) - h 1 ^ 10 ^) < 2 x 10" 47 , 

where the functions R r m {x), m > 10, are defined similarly to the ones given in the 
above theorem. We remark that in all the cases the maximum error is at x = 1. 
It is also important to notice that for bigger values of r we need bigger values of 
n in R r n (x) and h r n {x) to obtain similar bounds. 

Theorem 11.21 is simply a reformulation of the above theorem. 

4. On the time-parametrization of t c for the 
flsher-kolmogorov case 

By using normal forms theory it is well-known that in a neighborhood of the 
origin of system (jSJ) its unstable manifold can be parametrized as F(e rt ) for some 
analytic function F. So it seems natural to find bounds of the actual heteroclinic 
orbit which are rational functions of e rt . As far as we know this idea is new. 

First, we consider the family of rational functions 

= /9e rt + ae 2rt 

U 1 + (a + 2(3 - l)e rt + ae 2rt 

to try to approximate the function x r (t), where we denote by (x r (t) , y r (t)) the 
time-parametrization of the heteroclinic orbit r r . Notice that the parameters are 
taken in such a way that 

lim X(t) = 0, X(0) = -, lim X(t) = 1, 

t— oo 2 t— >oo 

properties that are also satisfied by x r (t). We do not care about the y component 
because y(t) = —x'{t). 

Recall that the vector field associated to system (jSJ) is 

Z r (x, y) = (Z[(x, y), Z r 2 (x, y)) = (-y, (r - l/r)y + x(x - 1)). 

To study the behavior of the flow of this system on the curve (X(t), Y(t)), where 
Y(t) = —X'(t) we compute 

M(t) := -Y'(t)Z[(X(t), Y(t)) + X'(t)Z r 2 (X(t), Y(t)). (10) 

If we introduce the compact notation $ := e rt we obtain that 

$ 3 (a(a + /3- 1)$ 2 + 2a$ + 0) P 3 ($) 



M{t) 



(a$ 2 + (a + 2/3-l)$ + l) 5 
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where P3 is a polynomial of degree 3 with coefficients depending also polynomially 
on a and (3. 

In order to simplify the expression of M(t) we take a — 1 — /?. Then 

$ 3 (2(l-/3)$ + /3)P 3 ($) 
1 J ((l-/3)$ 2 + /3$ + l) 5 ' 

Finally, taking /3 = 2\/2 — 2, the numerator and the denominator of the above 
fraction have a common zero. Then it writes as 

M(t) = ^17 + 12^(1 -6r 2 )$ 3 



(1 + v / 2 + $) 7 
In short we have proved the following: 

Theorem 4.1. Let (x r (t) , y r (t)) be the time-parametrization of the heteroclinic 
orbit T r of system ([8]) that satisfies x r (0) = 1/2. Define 

(2 + 2V2 + e rt )e rt 



[1 + V2- 



e 



rt\2 



Then: 

(1) When r < l/y/E it holds that sgn(x r (t) -X r {t)) = - sgn(t). 
(ii) When r = 1/y/E it holds that x r (t) = X'\t). 
(Hi) When r > 1/V6 it holds that sgn(x r (t) — X r {t)) = sgn(t). 

Theorem 11.31 follows from the above result, simply using that u = 1 — x and 
c = 1/r — r. 

Theorem [375] can also be used to obtain an explicit bound for x r (t). Since hi(x) = 

T 

—rx H x 2 < Y( x )j where (x,7 r (x)) is the parametrization of T r , solving 

2r 2 + 1 
the Cauchy problem 

x ' = rx ~ ^Vl x2 ' X(0) = 2' 

we obtain a function 

2r 2 + 1 

U r (t) 



1 + (4r 2 + l)e~ rV 

such that sgn(x r (t) — U r (t)) = — sgn(t). The other bounds given in Theorem 13.51 
give rise to implicit bounds of the form H(t,x) = of the curve (t, x r (t)). 

4.1. Sharper bounds when r £ (0, l/\/6). This subsection improves the re- 
sults of the previous section when r < 1/ \^6. 

Fixed r e (0, 1/V6) and given any natural n, following similar techniques that 
the ones used to prove Lemma I3~T| we can compute a function x r n {t) = X^=i ■> 
with $ = e rt and a\ = 1, that coincides with the solution x r (t) until order 2n in 
$. For instance 

02 = -77-5- — 7 and a 3 



2r 2 + l d (2r 2 + l)(3r 2 + 1)' 
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Looking $ as an independent variable, we compute the associated Pade approxi- 
mant of X]f=i a j^ > ' 7 °f order (n,n), obtaining 



n 

Z r n (t) = 2%*) = 3 4 , 

3=0 

where bj(r) and Cj(r) are polynomials on r. For instance 

3(2r 2 + l)(3r 2 + l)(4r 2 + 1)$ - 2(r - l)(r + l)(3r 2 + 1)$ 2 



3(2r 2 + l)(3r 2 + l)(4r 2 + 1) + 5(2r 2 + l)(3r 2 + 1)$ + (3r 2 + 2)$ 2 ' 



Notice that it satisfies that lim Z r n (f) = but we do not impose neither that 

t— >— oo 

Z;(0) = 1/2 nor that lim Z r n {t) = 1. 

t— >oo 



As in the proof of Theorem 14.11 to study the behavior of the flow associated 
to system (jHJ) on the curve (Z^(t), — Z^(t)') we compute the corresponding func- 
tion ffTOl) . It writes as 

M (t) = M ($) = - — 

n[ ) n{ ' (Qui®)) 5 

where Pk and Qk are polynomials of degree k with coefficients depending also 
polynomially on r. It can be seen that all these coefficients, as functions of r, 
take positive values on (0, l/y/E). We prove this fact introducing a new variable 
z satisfying r = z 2 /(v / 6(l + z 2 )) in each of the coefficients and then applying the 
Descartes' Theorem to the resulting polynomials. Then M n (t) > for all t and 
for all r in (0, 1/y/o). Moreover, this also shows that Z T n (t) is well defined because 

n 

Q n ($) = £ Cj (r)& > 0. 

3=0 

Since we want that Z r n (0) = 1/2 we need to modify Z^(t). To do this, for p > 0, 
we define the new family of functions W^(t, p) = Z^(p<&). It is not difficult to 
see that following the above procedure the same results hold. Moreover it can be 
proved that there exists a unique value po(r) G (0,1) such that W^(0, po(r)) = 
1/2. So we define X r n {t) = W^{t,p {r)) and it satisfies X*{0) = 1/2. Since the 
corresponding function ffTOl) is positive, for each n — 2, . . . , 8, it holds that 



sgn(x r (t)-X r (t)) = -sgn(t). 
Moreover lim X*{t) = b n (r)/c n (r) and it can be seen that these limits, which are 

t— >oo 

rational functions in r, satisfy 

1 < 6 8 (r)/c 8 (r) < 6 7 (r)/c 7 (r) < . . . < b 2 (r)/c 2 (r). 

Furthermore 

E k := max -—- - 1 

re(0,l/V6) V c fc( r ) 

decreases with fc. For instance E 2 = 3 — 4y / 5/3 ~ 0.02 and E 8 w 0.0007. 
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Conclusions 

In spite of the great interest of studying the traveling wave solutions of reaction- 
diffusion equations, u t = u xx + f(u), there are no methods for obtaining explicit 
bounds for them. In this paper we present an approach that allows to obtain 
this type of bounds in the general case. Moreover we introduce more elaborated 
tools, based on the control of the heteroclinic trajectories of an associated planar 
ordinary differential equation, that allow to improve the general results when we 
deal with a particular function /. We study with detail the classical Fisher- 
Kolmogorov case f(u) = u(l — u). The methods developed can be easily adapted 
to treat the Newell- Whitehead- Segel and the Zeldovich equations 

References 

[1] M. J. Ablowitz and A. Zcppetella, Explicit solutions of Fisher's equation for a special wave 
speed, Bull. Math. Biol. 41 (1979), 835-840. 

[2] R. A. Fisher, The wave of advance of advantageous genes, Ann. Eugenics 7 (1937), 355-369. 

[3] A. Gasull, H. Giacomini and J. Torregrosa, Some results on homoclinic and heteroclinic 
connections in planar systems, Nonlinearity 23 (2010), 2977-3001. 

[4] P. Grindrod, Patterns and Waves: The Theory and Applications of Reaction-Diffusion 
Equations, Clarendon Press (1991). 

[5] A. Kolmogorov, I. Petrovskii and N. Piscounov, A study of the diffusion equation with in- 
crease in the amount of substance, and its application to a biological problem, In V. M. 
Tikhomirov, editor, Selected Works of A. N. Kolmogorov I, 248-270. Kluwcr 1991. Trans- 
lated by V. M. Volosov from Bull. Moscow Univ., Math. Mcch. 1 (1937), 1-25. 

[6] A. C. Newell and J. A. Whitehead, Finite bandwidth, finite amplitude convection, J. Fluid 
Mech. 38 (1969), 279-303. 

[7] F. Sanchez-Garduho and P. K. Maini, Travelling wave phenomena in some degenerate 
reaction- diffusion equations, J. Differential Equations 117 (1995), 281-319. 

[8] F. Sanchez-Garduno and P. K. Maini, Travelling wave phenomena in non-linear diffusion 
degenerate Nagumo equations, J. Math. Biol. 35 (1997), 713-728. 

[9] L. A. Segel, Distant sidewalls cause slow amplitude modulation of cellular convection, J. 
Fluid Mech. 38 (1969), 203-224. 
[10] J. Xin, Front propagation in heterogeneous media, SIAM Rev. 42 (2000), 161-230. 
[11] Y. B. Zeldovich and D. A. Frank-Kamenetskii, A theory of thermal propagation of flame, 
Acta Physicochimica URSS 9 (1938), 341-350. English translation: In Dynamics of curved 
fronts, editor R. Pelce, Perspectives in Physics Series, Academic Press, New York (1988), 
131-140. 

Dept. de Matematiques. Universitat Autonoma de Barcelona. Edifici C. 08193 
Bellaterra, Barcelona. Spain 

E-mail address: gasullOmat . uab . cat 

Laboratoire de Mathematiques et Physique Theorique. Faculte des Sciences 
et Techniques. Universite de Tours. 37200 Tours. France 
E-mail address: Hector.Giacomini@lmpt.univ-tours.fr 

Dept. de Matematiques, Universitat Autonoma de Barcelona, Edifici C. 08193 
Bellaterra, Barcelona. Spain 
E-mail address: torre@mat .uab . cat 



